Electron turbulence at nanoscale junctions 
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Electron transport through a nanostructure can be characterized in part using concepts from 
classical fluid dynamics. It is thus natural to ask how far the analogy can be taken, and whether 
the electron liquid can exhibit nonlinear dynamical effects such as turbulence. Here we present 
an ab-initio study of the electron dynamics in nanojunctions which reveals that the latter indeed 
exhibits behavior quite similar to that of a classical fluid. In particular, we find that a transition 
from laminar to turbulent flow occurs with increasing current, corresponding to increasing Reynolds 
numbers. These results reveal unexpected features of electron dynamics and shed new light on our 
understanding of transport properties of nanoscale systems. 
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Electron transport through a nanoscale junction is usu- 
ally described as a scattering problem i 1 ' 2 ' 3 ' ' 5 ' 6 ' On the 
other hand, it has been shown^ i 8 ' 9 ' 10 that the behavior of 
the electron liquid obeys dynamical equations of motion 
which arc similar in form to those governing the dynamics 
of classical liquids. In particular, it was recently shown 10 
that the time-dependent Schrodinger equation (TDSE) 
for electrons flowing across a nanostructure can be cast 
in the form of generalized Navier-Stokes equations. A 
consequence of this analogy is the prediction that under 
certain conditions the laminar flow of the liquid may be- 
come unstable and turbulent behavior is expected i 10 ' 11 ! 12 
One can then borrow knowledge from classical fluid dy- 
namics and hypothesize that the electron flow will make 
a transition from laminar to turbulent regimes, if, e.g., 
the current is increased. Like in a classical liquid, if, for 
instance, electrons flow from one electrode (call it the 
top electrode) to another (call it the bottom electrode) 
across a nanojunction, turbulence would first manifest 
itself with the break-up of the top-down symmetry of the 
electron flowi 11 ' 12 By increasing the current further one 
should observe the formation of eddies in proximity to the 
junction. These eddies are created because of the larger 
kinetic energy in the direction of current flow (longitudi- 
nal kinetic energy) compared to the transverse direction 
(transverse kinetic energy). By increasing the current 
further, the disparity between longitudinal and trans- 
verse components of the kinetic energy increases, and 
the system flow eventually breaks any remaining sym- 
metry, thus fully developing turbulencei 11 ' 12 Despite the 
prediction of turbulent behavior of the electron liquid in 
nanostructuresi - an explicit demonstration of this phe- 
nomenon and the analysis of its microscopic features has 
not been presented yet. 

In this communication we set to show numerically tur- 
bulent effects for the electron liquid crossing a nanojunc- 
tion both by solving directly the TDSE, and by solv- 
ing the generalized Navier-Stokes equations derived in 
Ref. G1. These read 



An(r,f) 
mn(r,t)D t Ui(r,t) 



-n(r,i)V • u(r, t) , 
-diP(r,t) + djirij 
-n(r, t)diV ext (r,t) . 



Here, n is the electron density, u = j/n is the velocity 
field, i.e. the ratio between the current density j and the 
density, and D t = Jj+u- V is the convective derivative.— 
P(r, t) is the pressure of the liquid, I4 x t(r, t) an external 
potential, and m,j is a traceless tensor that describes the 
shear effect on the liquid. It has the form 



(2) 



7Tj,j = // ( diUj + djUi - -5ijd k Uk 



where rj — hnf(n) is the viscosity of the electron liquid 
and f(n) is a smooth function of the density^ 

In analogy with the classical case we expect that the 
atomic structure, and in particular atomic defects in 
proximity to the junction, play the role of "obstacles" 
for the liquid and thus favor turbulence. Since we aim at 
showing that turbulence develops irrespective of the un- 
derlying atomic structure, we consider electrons interact- 
ing with a uniform positive background charge (i.e., the 
"jcllium" model 14 ). The system we consider therefore 
consists of two large but finite jellium electrodes- sub- 
ject to a bias- connected via a nanoscale jellium bridge. 
(The jellium edge of this system is represented with solid 
lines in each panel of Fig. [TJ) We choose the density of 
the jellium at equilibrium typical of bulk gold (r s ~ 3ao). 
For computational convenience we choose a quasi-2D sys- 
tem, approximately 2.8 A thick^ Note that, everything 
else being equal, a quasi-2D geometry disfavors turbu- 
lence compared to a 3D one. We thus expect that if tur- 
bulence develops in our chosen quasi-2D geometry with 
given thickness, then turbulence will develop even more 
easily if we leave everything else unchanged (including 
the total current) and increase the thickness of the elec- 
trodes. 

The solution of the TDSE for the many-body system is 
obtained within Time-Dependent Current Density Func- 
tional Theory (TDCDFT); 8 i.e., for each single-particle 
state <p a , we have solved the equation of motion (in 
atomic units) 
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where c is the speed of light, Vj e \ is the potential due 
to the jellium, i>h is the Hartree potential, and v xc is the 
exchange-correlation scalar potential.— The shear viscos- 
ity of the electron liquid r/ enters the problem through the 
exchange-correlation vector potential a xc j&i£ If we make 
the appoximation that, at any given time, the viscosity 
is a function of the density only, and does not depend 
on time explicitly, we find that the i th component of a xc 
evolves according toi^ 

c dt ~ n^Y 1 dr 2 3 did, 
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For the viscosity r\ of the electron liquid we have used 
the one reported in R ef. IT3. We employ the approach 
described in Rcfs. to initiate electron dynam- 

ics and calculate the current. We prepare the system by 
placing it in its ground state; because of the applied bias, 
the system exhibits a separation of charge. At a time 
t — 0, we remove the bias, and let the system evolve ac- 
cording to eq. 31 After long time scales, the electrons 
encounter the far boundaries of the jellium electrodes 
and reflect. However, we are interested in comparing 
the electron dynamics calculated from eq. {3]) with the 
one obtained by solving the Navier-Stokes equations fTJ) , 
at times smaller than this reflection time. Equations ((T|) 
are solved by assuming the same density, initial veloc- 
ity and viscosity of the liquid employed in the TDCDFT 
calculation^ We also solve eqs. ((TJ assuming the liq- 
uid to be incompressible. This simplifies the calculations 
enormously but leads to some minor differences with the 
solutions of eq. §S§ (see discussion later). 

Fig. [1] depicts the flow of electrons across the nanos- 
tructure, for a range of biases between 0.02V and 3.0V, 
after the initial transient. The panels (a)-(d) correspond 
to the solution of eq. ©; the panels (e)-(h) to the solution 
of the eqs. (fTJ) using the same set of parameters^ Panel 
(a) has to be compared with panel (e); panel (b) with (f), 
and so on. As anticipated, we observe some differences 
between the solutions of the eqs. |T]) and the solutions of 
the eq. Q. These differences are due to the details of the 
charge configuration at the electrode-junction interface, 
and some degree of compressibility of the quantum liquid 
in the junction. 

From Fig. Q] we can see the effect of surface charges, 
in that some electrons flow parallel to the surfaces^ 
More importantly, at low biases, the flow is laminar and 
"smooth" . In addition, at these biases the current den- 
sity shows an almost perfect top-bottom symmetry: the 
direction of the flow is symmetric with respect to the op- 
eration z — > — z. This symmetry is even more evident by 
comparing the curl of the current density in the top and 
bottom electrodes (see for instance Fig.QJa) and (e)). 

By increasing the bias, however, a transition occurs: 
the symmetry z — > — z of the current density breaks com- 
pletely, and eddies start to appear in proximity to the 



junction. This is clearly evident, for instance, in Fig.QJd) 
and (h). The outgoing current density in the bottom elec- 
trode has a more varied angular behavior, in contrast to 
the behavior in the top electrode, in which the electron 
liquid flows more uniformly toward the junction. 

Since the panels (e)-(h) of Fig. Q] practically describe 
the dynamics of a classical fluid with the same parameters 
as the quantum liquid, the analogy between the electron 
flow and the one of a classical liquid is quite evident. We 
can push this analogy even further by defining a Reynolds 
number for the quantum system as well: R — u z Lp/rj, 
where u z is the longitudinal velocity in the junction, L 
is the width of the junction, and p is the density. Using 
the density of valence electrons in gold, and using the 
current density in the junction at t = 1.4 fs, we obtain 
for the quantum case the following Reynolds numbers: 
0.216, 2.16, 10.8 and 32.5 for 0.02 V, 0.2 V, 1.0 V, and 
3.0 V, respectively. 

Just like in the classical case, we can then reinterpret 
the above results as follows. At low Reynolds numbers, 
the flow is highly symmetric from top to bottom. This 
symmetry is lost as the Reynolds number is increased. 
At high Reynolds numbers, the incident flow is lami- 
nar, while the outgoing flow has a jet-like character, and 
"turns back" on itself creating local eddies in the current 
density. 

Having shown the similarity between the current flow 
obtained using the Navier-Stokes equations ([I]) and the 
one obtained solving eq. we can study the first one at 
times scales prohibitive for full quantum mechanical sim- 
ulations. We can also study the effect of a larger thick- 
ness of the electrodes on the turbulent flow by realizing 
that this is equivalent to increasing the Reynolds num- 
ber (which, incidentally, is also equivalent to increasing 
the bias). This is illustrated in Fig. [5] where the current 
density and the curl of the current density are plotted 
for the Reynolds number 32.5 of Fig. [TJh) (left panel of 
Fig. [5]); same system but with a Reynolds number five 
times larger (middle panel of Fig.fJ]); and (right panel of 
Fig. [2]) with a Reynolds number ten times larger. 

From Fig. [5] it is evident that by increasing thickness 
the last remaining symmetry x — ► — x is broken at ear- 
lier times, leading to turbulent behavior closer to the 
junction. For instance, in the case represented in Fig. [2] 
(middle panel), the left-right symmetry is lost at about 
14 fs, with consequent asymmetric flow within about 50 
A from the junction center. For the structure repre- 
sented in Fig. |2] (right panel), the symmetry is broken at 
about 6 fs, and the flow asymmetry appears at about 25 
A from the junction. 

We can better quantify the amount of turbulence by 
calculating the velocity correlation tensor— 

B ik = ((«j(r) - Vi(r + Sr))(v k (r) - v k {v + Sr))) , (5) 

where Sr is a given distance, and i, k = x, y, z. Here, the 
angle brackets (• ■ •} denote averaging over all positions r 
within a given region. Fully developed and isotropic tur- 
bulence has a velocity correlation tensor that is a function 
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FIG. 1: (Color online) Panels (a)-(d): Electron current density for electrons moving from the top electrode to the bottom 
electrode across a nanojunction at t = 1.4 fs, for an initial bias of (a) 0.02 V, (b) 0.2 V, (c) 1.0 V and (d) 3.0 V. The arrows 
denote the current density, while the level sets denote the curl of the 2D current density. The solid lines delimit the contour of 
the junction. Panels (e)-(h): Velocity field solution of the equations ([T]), for a liquid with same velocity, density and viscosity 
as the quantum mechanical one. 
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FIG. 2: (Color online) Current density (arrows) and curl of 
the current density (denoted by level sets) of the electron 
liquid, for three different Reynolds numbers, 32.5 (left panel), 
162 (middle panel), and 325 (right panel). Note that the fluid 
velocity has lost perfect left-right symmetry in the middle- 
and right-panel cases. 



only of the magnitude of <5r, and increases quadratically 
with distanced Instead, the turbulence in the examples 
of Fig. [J] is not fully developed. The velocity correlation 



tensor, Bik, thus depends on both the magnitude of 5r 
as well as its direction. This is illustrated in Fig. [3] where 
various components of Bik are plotted at t = 75.0 fs for 
the system with Reynolds number 32.5 as a function of 
the magnitude of 6r, where we have chosen Sr to point 
in the longitudinal (z) direction. The spatial averaging 
has been carried out over the left-hand side of the out- 
going region (that is, in the region z = [-27.3 A, -72.1 
A],x= [-25.9 A, 0.0 A], where the origin is in the center 
of the junction). For comparison, the same quantity is 
plotted for the laminar case, i.e. for a Reynolds number 
of .216. (To compare the laminar and turbulent cases we 
have scaled the average turbulent velocity to the average 
laminar velocity.) As expected, in the laminar case the 
correlation tensor is essentially zero, while for the turbu- 
lent case it increases with distance. 

We conclude by noting that an experiment in which 
the electron flow can be monitored directly may measure 
nonlaminar electron behavior as an asymmetry between 
the incoming and outgoing patterns of the current den- 
sity through a nanojunction. Experiments similar to the 
ones reported in Ref. [20l which use scanning probe mi- 
croscopy to image the flowlines, may provide such capa- 
bilities. Note, however, that since these scanning probe 



4 



0.2 
0.15 
"J 0.1 

°< 



1 


i i i 


















Laminar 

i i i i 



3 6 9 12 15 

d [A] 

FIG. 3: (Color online) Various components of the velocity 
correlation tensor Bit as a function of distance d = | <5r | , for a 
Reynolds number of R = 32.5. For the case where R = 0.216, 
the elements of Bit are orders of magnitude smaller, and so 
the corresponding curves for the laminar case coincide with 
the x axis. 

techniques record images over time scales much longer 
than the turbulent electron dynamics, one would expect 



that in the turbulent regime images of current flow ap- 
pear "smeared out" compared to the laminar case, and 
asymmetric in the incoming and outgoing patterns. From 
the pre sent work and the analytical results obtained in 
Ref. [l(J, we also suggest that fully 3D and non-adiabatic 
junctions, i.e. junctions with a geometry that changes 
abruptly (like the one explored in this work), are the 
best candidates to observe turbulence. We also expect 
defects and other impurities to favor turbulent behavior 
by playing the role of "obstacles" for the electron flow. 

We finally note that turbulent behavior may have 
consequences on the formation (or lack thereof) of 
local equilibrium distributions in the electrodes^ and 
may generate non-trivial local electron heating effects 
in the junctional and at the eddies sites; phenomena 
and properties which are still poorly understood at 
the nanoscale and ultimately may have unexpected 
consequences on the stability of nanostructures under 
current flow.— 
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